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Abstract 

We study the peculiarities of the nonstationary Casimir effect (creation of photons in cavities 
with moving boundaries) in the special case of two resonantly coupled modes with frequencies loq 
and (3 + A)ui , parametrically excited due to small amplitude oscillations of the ideal cavity wall 
at the frequency 2u>o(l + S) (with \S\, |A| <§; 1). The effects of thermally induced oscillations 
in time dependences of the mean numbers of created photons and the exchange of quantum 
purities between the modes are discovered. Squeezing and photon distributions in each modes 
are calculated for initial vacuum and thermal states. A possibility of compensation of detunings 
is shown. 
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1 Introduction 



Classical and quantum phenomena in cavities with moving boundaries attracted attention of many 
researchers for a long time (see review Q). Especially popular this topic became in the last 
decade, being known now under the names nonstationary Casimir effect dynamical Casimir 
effect ||], or mirror (motion) induced radiation ||]. One of several theoretical results obtained 
in the last years was the prediction of the exponential growth of the energy of the field under the 
resonance conditions, when the wall performs vibrations at the frequency which is a multiple of the 
unperturbed field eigenfrequency (B|-|8| . 

In most papers, the special case of the one-dimensional cavity with an equidistant spectrum 
of unperturbed field modes was studied. Then an infinite number of modes are excited due to an 
intermode interaction (resulting from the Doppler effect on the moving boundary), and the total 
number of created photons depends on time as t 2 (but the total energy increases exponentially) 
]l], |8|-|lT|] . A more realistic case of a three-dimensional cavity was considered, e.g., in |T2]-|l4|. In 
particular, simple analytical solutions describing the resonance creation of photons in cavities with 
totally nonequidistant spectra, when only one mode is in resonance with the vibrating wall, were 
given in [|],[7|,^]. However, the case of a few resonantly interacting modes is also of a great interest. 
In particular, the interaction between the excited field mode and the detector, approximated either 



by a harmonic oscillator or by a two-level atom, was considered in |7],p, 15-1 



It was shown recently [20|, that for rectangular three- and two-dimensional cavities, there exist 
special configurations, when two (or more) modes can occur in resonance with the moving wall 
and between themselves. For example, the eigenmode spectrum in a cubical cavity of length L is 
given by the formula u>kim = (ttc/L)\/ k 2 + I 2 + m 2 . Since W511 = 3cjiu, both the modes, {111} and 
{511}, will be excited if the wall oscillates at the frequency 2wm. In such a case, it appears |20| that 
the energies of both the modes increase in time exponentially (with some additional oscillations), 
but the increment of the exponential growth is twice smaller than in the case of a single resonance 
mode. 

Our goal is to give analytical solutions describing the case of two resonantly interacting modes 
in a more or less generic situation, without specializing the form of the cavity. We take into 
account the possibility of detuning between the frequency of the vibrating wall and the twice 
fundamental field eigenfrequency, and detuning between the frequency of the second mode and the 
triple fundamental eigenfrequency (however, we do not consider the effects of damping). Using 
the method of slowly varying amplitudes, we obtain explicit solutions of equations of motion for 
the canonical operators describing the fields in coupled modes, which enable us to calculate the 
mean energy and the invariant uncertainty product for each mode. We analyze the evolution of 
the photon distribution functions for each mode and show the possibility of strong squeezing in the 
long-time limit (for other configurations this was done in |7|,|8|,|l(|). 
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2 General solution for two coupled modes 

We use the Hamiltonian approach proposed by Law (|| and developed in |2lJ (for other references 
see H]). Consider the scalar massless field <3?(r, t), satisfying the wave equation = V 2 <5 inside 
the cavity and the Dirichlet boundary condition = on the boundary. We assume that we know 
the complete orthonormalized set of eigenfunctions (and eigenfrequencies) of the Laplace equation 
V 2 / a (r) + unfair) = in the case of stationary cavity. Now suppose that a part of the boundary 
is a plane surface moving according to a prescribed law of motion L(t) (for the most recent study of 
the case when L(t) is a dynamical variable due to the back reaction of the field see (2^]). Expanding 
the field <&(r,t) over "instantaneous" eigenfunctions f a (r;L(t)), 

*(r,t) = 5>a(*)/a(r; £(*)), (1) 

a 

we satisfy automatically the boundary conditions. Then the dynamics of the field is described 
completely by the dynamics of the generalized coordinates q a (t), which, in turn, can be derived 



from the time- dependent Hamiltonian [21| 



= \ E [pI + <&W))&[ + T77Y E V«m am (2) 
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with antisymmetrical time-independent coefficients 



m al3 = -mp a = L J dV 9 ^ a J^ L fp(r;L). (3) 

For example, in the case of a rectangular three-dimensional cavity, the eigenmodes are the well 
known products of three sine functions like sin (irk x x/L x ), labeled by three natural numbers 
k x ,k y ,k z . If one surface of the parallelepiped, perpendicular to the x-axis, moves in the x-direction 
(so that the L^-dimension of the cavity is a function of time) , then |2(J 

m ki = (-l) k ^^S kyjy S k ^. (4) 

Jx ™x 

We are interested in the case when one of the cavity walls performs small oscillations with the 
frequency Q, close to the double frequency of some unperturbed mode ui^ = 1 (i.e., we normalize 
all frequencies by u\ ), so that the time-dependent frequency 0J\{t) reads 

ui(t) = 1 + 2ecos(2cJt), uj = 1 + 5, (5) 

where we assume that \5\ <C 1 and |e| <C 1. Also we suppose that the unperturbed field frequency 
spectrum includes the frequency cdj = 3 + A with |A| <C 1, but it does not contain frequencies 
close to 5o»j . Then we have the case of two resonantly interacting modes, and it is sufficient to 
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consider only the part of the total Hamiltonian (0) related to these modes: 



flia = -(p? + p§) + -[l + 4ecos(2wt)]a:? 

+-[9 + 6A + ecos(257t)]:p§ 
+3/ie sin(2u;i) (pix 3 - PzX\) . 



(6) 



The constant parameter [i is proportional to the coefficient m\2 in (|2|). Writing (|6|) we have 
neglected the second order terms with respect to e and A. Since the time-dependent frequency 
shift u>i(t) — u>i ~ L(i) — Lo is chosen to be proportional to cos(2uJt), the last term in (||), which is 
proportional to L(t), must depend on time as sin(2uJt). All numerical coefficients in @ are chosen 
in order to avoid an appearance of fractions in formulas below. Parameter e has the same order of 
magnitude as e, but it does not affect the solution in the zeroth order approximation (in which we 
are interested here), as will be shown below. 

Hamiltonian (|6j) results in the following differential equations for the generalized coordinates x\ 
and X3 (we neglect corrections of the second order): 



x\ = — [1 + 4e cos(2u;t)] x\ 

+24//e [cos(2tJt)x3 + sin(2u;t)i;3] , 



(7) 



x 3 = - [9 + 6 A + e cos(2uT£)] x 3 

—2A/ie [cos(2u;t)xi + sin(2aJt)i;i] . 

We solve equations (|7|) and (||) , using the method of slowly varying amplitudes 
applied earlier in studies 0-^,^]). Namely, we look for the solutions in the form 



x k (t) = $(t)e ik " t + ^(t)e 



-ikujt 



k= 1,3 



(8) 

(which was 
(9) 



where each coefficient £1. is a slowly varying function of time, whose derivatives are proportional to 
the small parameters e, 5, A, so that we can neglect the second-order derivatives £ fc . Then, equating 
the coefficients at exp[±iuJ£] in equation ([?]) and coefficients at exp[±3iuJt] in equation (|J) (and 
neglecting the terms proportional to squares and products of small coefficients e,5, A), we obtain 
a set of equations with constant coefficients for the slowly varying amplitudes. It is convenient to 
write this set in the matrix form dw/dt = Av, introducing vector v = ^£1% £i~> £3") an d matrix 



A 



-18 


ie 


12ifie 





—ie 


id 





— 12ifie 







i(A - 35) 








— 4i//e 





-t*(A - 35) 



(10) 



Strictly speaking, function x\{t) contains also terms proportional to exp[±3iu>£], as well as function 
x 3 (t) contains terms proportional to exp[±iuJt]. However, the amplitudes of these additional terms 
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are proportional to the small parameters e, 5, A, so we neglect them, confining ourselves to the 
principal terms whose amplitudes are of the zeroth order with respect to small parameters. Just 
for this reason, the parameter e does not appear in matrix A: while the term 4ecos(2tJt) standing 
at x\ in equation ([?]) connects the slowly varying amplitudes £± and a similar term ecos(2tJi) 
standing at x% in (||) does not connect the amplitudes £3 and £3" between themselves, but it gives 
small corrections to the amplitudes of higher-order terms in X3 , oscillating with frequencies U and 
5tU. A more strict (from the mathematical point of view) method, based on the multiple scale 



analysis, was exposed in [20]. In the zeroth-order approximation it leads to the same results. 
Matrix A has four eigenvalues, ±A±, where 



1 

1 
2 



6 2 (1 



a ± \fc 
+ Vb± 



Vb 



5 2 - (A-35) 2 , 



(11) 
(12) 

(13) 



= 2(<S-e)(A-3<5) + ve 2 
x 2(<J + e)(A-3«5) + ue 2 



(14) 



c = a 



+ 



2e 2 v 



(A - \S) 2 - e 5 
2 



e z + {A - 26){A - 46) 



96fi 2 



(15) 
(16) 



Comparing our matrix A ( |10|) with a similar matrix found in [gOJ for the rectangular cavity, one can 
verify that if modes {k x ,m, n} and {j x ,m, n} are in resonance, then [i = j x / (12k x ) (different phases 
of elements of matrices in and here are due to different choices of trigonometrical dependences 
in function L(t): we use cos- function instead of sin- function in pCf ). In particular, for the modes 
{111} and {511} of the cubical cavity we have v = 50/3. Due to this explicit example (and other 
examples related to the rectangular cavities), we assume hereafter that parameter v is large, so 
that it satisfies the inequality v 3> 1. 



3 Exact resonance 



We start with the simplest case of the exact resonance, S = A = 0. Then formula (0) yields 

a£> = I (1 ± VT^) . (17) 

If v = (no intermode coupling), then we obtain the increment of the exponential growth = e, 
in accordance with the solutions of the single-mode problem found in [@, | , 26 1 . On the contrary, 
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for v > 1/2 the increment (real part of A±) does not depend on the form of the cavity (which is 
"hidden" in the value of i>), being exactly twice smaller than in the single-mode case. In the special 



case of rectangular cavities this result was obtained in [20]. 

After some straightforward calculations we arrive at the following explicit expressions for the 
time dependences of the generalized coordinates (as far as we neglect all corrections of the order 
of e in the amplitude coefficients, we may identify the canonical momenta with velocities, = £&, 
k = 1,2): 

sin(pr)" 



xi (t) =xi(0) 



C 1 cos(/ot) + S 1 



p 



-Pi(0) 



S 1 cos(pr) + C l 



sin(pr) 



+8/i 



sm(pr) 
P 



3Srs 3 (0) + CrP3(0) 



(18) 



where 



x 3 (t) = x 3 (0) 



(7 3 + cos(pr) - S 3 



sm(pr) 



P 



St cos(pr) - C 3 



sin(pr) 



-8/i 



sin(pr) ■ 
P 



S+ Xl (0)-C+ Pl (0) 



t = -et, p = \J1v- 1, 



(19) 



(20) 



C^(r;t) = cosh r cos (ku> t) ± sinhr sin(/cu;t), 
S^(r;t) = sinhr cos (k&t) ± cosh r sin (kut). 

Similar expressions for the time dependent canonical momenta can be easily obtained by differen- 
tiating equations (|lS|)-(|l9|) with respect to the "fast time" t, considering the "slow time" variable 
r as an independent parameter and neglecting corrections of the order of e, 5, A in the amplitude 
coefficients. The following relations hold in this approximation: 



9Ct 



±kS^, 



9S£ 



at k ' at 

Symbols Xk,Pk in equations (|l8|)-(19) can be considered both as classical variables and quantum 
operators in the Heisenberg picture, due to the linearity of the problem (or due to the quadratic 



nature of Hamiltonian (§)). Using equations dl8|)-(19) and their momentum counterparts, one can 
calculate mean values of squares and products of canonical variables (operators) at any moment 
of time, provided such mean values were known at the initial moment t = 0. We confine ourselves 
to the simplest case, when initially the field modes were in thermal states with the mean photon 



6 



numbers (0i — l)/2 and (#3 — l)/2, where 9k = coth(A;/3/2), (3 being inverse absolute temperature 
in dimensionless units. One can check the relations 



731 = ^ 



r 1 

; i3 - ^2 



The mean energies in each mode, = ^(pf + u)f.xj,), depend on time as follows 



f + 3 



-\ cosh(2r) 



sin 2 (pr) 



+ cos 2 (pr) + sinh(2r 



sin(2/3r) 1 



(21) 



fa 



30, 



cosh(2r) 



sin 2 (/)r) 



(1 + 21/013) 



+ cos 2 (/>r) — sinh(2r 



sin(2pr) 1 



(22) 



P J 

For rectangular cavities (and, perhaps, for others), v 3> 1 in the case of intermode resonance. Then 



P 



2v 3> 1, so that equations (f2~I|) and ( |22"| ) can be simplified: 



fi 
£3 



1 

2 
3 



cosh(2r) 0i cos 2 (/or) + 63 sin 2 (pr) 
cos1i(2t) 03 cos 2 (/or) + 0i sin 2 (pr) 



For the initial vacuum states 
of photons in each mode: £k ? 



?1 = 03 = 1) we have a monotonous growth of energy and number 
(a»l°^/2) cosh(2r), k = 1,2. On the contrary, for high-temperature 
initial states (with equal temperatures), 0i = 303 3> 1, which results in strong oscillations of mean 
energies and numbers of photons: 

2 



£1 
£3 



— cos1i(2t) 



1 



■ sin (pr) 



yCosh(2-r) l + 2sin 2 (pr) 



In Fig. 1 we show normalized time dependences of the mean energies (calculated with the aid of 
complete formulas ( ^l| ) and (^)) in each resonant mode in the low- and high-temperature limits. 
Note that one has 0i « 140, if L = 1 cm and T = 300 K. 

Since Hamiltonian (|2|) is quadratic with respect to canonical coordinates and momenta), the 
initial thermal state is transformed with time to a generic Gaussian quantum state (whose density 
matrix or Wigner function is a Gaussian exponential). The single-mode density matrices of such 
states are completely characterized (in the case of the mean values of quadrature components), 
besides the mean energy, by the invariant uncertainty product (IUP) 



V = {x 2 )(p 2 ) - {(xp + px)/2f 



(23) 
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For v > 1/2 this quantity varies in time periodically. For the first mode, 

.21/031-1 



Vs 



cos (pr) + sin (2pr) 



2(2z/ - 1) 



+ sin 4 (/?r) 



21/031 + 1 



2i/- 1 

For another excited mode one should interchange indices 1 and 3 in (]2 



The 'purity of quantum Gaussian states is expressed through IUP as Tr(/> 2 ) 
is the statistical operator of the state). For initial vacuum states, 



V 



1,3 



1 + 



■ sin 4 (pr) 



(24) 



(4P)-V2 (here 



(25) 



(2i/ - l) 2 

so that the states remain practically pure, if v S> 1. This is a significant difference from the case 
of a single resonance mode coupled to a harmonic oscillator detector, studied in ]7],||]. If 9% ^ 03 
(nonzero initial temperature) and v S> 1, then 



01 cos 2 (/rr) + 03 sin 2 (pr) 



so that coupled modes exchange their purities at the moments when sm(pr) = 1. 

The squeezing coefficient, defined as the ratio of the minimal value of any canonical variance (o~ x 
or o~p) for the period of fast oscillations (with frequencies ZJ or 3uJ) to the vacuum value ((2^^ ) _1 
or w^/2), can be expressed through £ = £/&% and V as |27|, p8|| 



2V 



(26) 



£ + V£ 2 - V 

For thermal states, s(0) = 9, but asymptotically one can obtain any desired degree of squeezing 
(even for initial high-temperature states) in each resonant mode, because s ~ D(t)/£(t) for t>1. 
Evidently, £ is nothing but the mean number of photons in the mode, if £ 3> 1. The variance of 
the photon distribution in the Gaussian states is given by the formula [29| 



2£ 2 -V- 1/4. 



The photon statistics is strongly super-Poissonian for r> 1, when a n /{n) ~ 2£ ^> 1. Nonetheless, 
the quantum states of each excited mode become highly nonclassical for r> 1. It is seen from the 
photon distribution function (PDF), which can be expressed in terms of two parameters, £ and T>, 
for any Gaussian state with zero mean values of quadrature variables |29||): 



1 + 4V 



A£V /2 



VI + 4£ + 42? \l + 4V + 4£ i 
P • 42? - I 

J 97. 



; (42? + l) 2 - 16f 2 



(27) 
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Here P n (z) is the Legendre polynomial. Initially c (0) = X>(0), and (27) is a monotonous geometric 
distribution, V n (0) = 2(9 — l) n /(9 + l) n . Since we are interested in the cases of large numbers of 
photons created due to the parametric resonance, it is convenient to use asymptotical forms of the 
exact distribution (|27| ) for n> 1, Note that the argument of the Legendre polynomial in (|27| ) is 
always outside the interval (—1,1), being equal to 1 only at the initial moment (for thermal states). 
With the course of time this argument increases to oo, becomes pure imaginary, and asymptotically 
goes to zero, when £ 3> T>. Therefore it is convenient to use the asymptotical formula for n ^> 1 [[30| 

/ c \l/2 

P n (cosh£) « -4- J ( n + 1/2] £) (28) 



V sinh £ , 

(where Iq(z) is the modified Bessel function), because it holds even for complex values of variable £, 
provided Re£ > and |Im£| < ir. We are interested in the case when the mean energy of each mode 
significantly exceeds its initial value, i.e., £ S> y/V. If V 3> 1 (high temperature initial states), 
then it is possible that £ < T> + 1/4. In such a case, £ is real and large, so that one can replace the 
modified Bessel function Iq(x) by its asymptotical expression for Iq{x) (2-kx)~ 1 I 2 exp(x). 

Moreover, one can use either of approximate equalities exp(x) ~ 2cosh(x) ~ 2sinh(x). 

When £ > D + 1/4 (this is just the regime when squeezing happens), the argument of the 
Legendre polynomial becomes pure imaginary, and the variable £ becomes complex: £ = —m/2 + y, 
with y > 0. Then one should replace the modified Bessel function in ( p8[ ) by the usual Bessel 
function Jo ([n + 1/2] [tt/2 + iy]). Since the absolute value of the argument of this function is also 
large, one can replace Jq{x) by (irx/2)~ 1 / 2 cos(x — vr/4). In this case, we have different expressions 
for even and odd values of index n (this is especially clear from the initial formula (p7|), when the 
argument of the Legendre polynomial is close to zero for £ 3> T>). Finally, we arrive at the following 
asymptotical expression, which holds both for real and imaginary values of the argument of the 
Legendre polynomial in (|27|) (provided £ 2 3> T>): 

n/2 



V2 



l + 4V-A£ 



Vwn£ (l + 4£> + 4<S 

J cosh ([n + 1/2] In Ixl) , n even 
I sinh ([n + 1/2] In Ixl) , n odd 



(n+l)/2 

(29) 



4P-1 + W£ 2 - V 
X = ; . • (30) 



where 

47 

X = — 

' \AV + l) 2 - \Q£ 2 

If the ratio £/T> is not very large, then In \x\ is not small, and one can replace cosh and sinh functions 
in ( p9|) by exponentials. In this case we have a nonoscillating "quasigeometric" distribution, i.e., a 
geometric distribution modified by a slowly decreasing factor n -1 / 2 : 

n+l/2 

(31) 



1 UV-l + 4\/£ 2 -V\ ' 

n ~ T^J v 4p + 1 + 4f ; 
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On the contrary, if £ S> T> (then \ is pure imaginary and small), we have a strongly oscillating 
distribution for n < £ /V (typical for squeezed states): 

, 4£> + l\ n / 2 

V n ~ \ ~ 1 



7rn£ V 2£ 

cosh ( n 4 ^ 1 ) , n even 
smn In ^ I n odd 

Only the "tail" of distribution (|32|) does not oscillate (and does not depend on V): 



exp 

Vn « 



-n/(2f) 



425-1 

n ^- > 1. (33) 



V2vrn? 45 

For the vacuum initial state, the PDF oscillates from the beginning. For r>l and i/> 1 we can 
write (for each excited mode) 

n/2 



V n « e -J-L(l-4e-*V 
V 7rn V / 

cosh sin 4 (/9r)e _2r ^ , n even 
sinh ^ sin 4 ( / or)e~ 2 ' r J , n odd 

Oscillations exist for any n at the moments of time when sin(pr) = 0. Two examples of the photon 
distribution functions for the initial vacuum and thermal states (calculated using the exact formula 
(|27|)) are given in figures 2 and 3. 



4 Influence of detunings 

Now let us consider the generic case of nonzero detunings. Hereafter we use the normalized pa- 
rameters 5 = 5/e and A = A/e. The condition of the parametric resonance is that at least one of 



the eigenvalues A±, given by equations (11) or (|12|), has nonzero real part. Since we assume that 
v > 1, the coefficient a (13) is always negative. Looking at the eigenvalues in the form (|i"l|), we see 
that there are two possibilities. 

The first one is to choose c < 0. Then complex number a + yfc has nonzero imaginary part, 
so its square root inevitably has nonzero real part. Designating A — 4<5 = r/, we see that c < for 
\r}\ < rjc, where the critical value ry c (which depends on v and 5) must be obviously less than 1. If 
v S> 1 and \S\ ~ 1, then r/ c is close to 1, with corrections of the order of v~ l . If |<5| 3> 1, then for 
1 ~ \r]\ <C |£| one can write c ~ 2v {rf — 1) +4<5 2 n 2 (taking into account that A — 25 = 25 + r] ~ 25). 
In this way we obtain the following approximate inequality, describing the region of resonance in 
the parameter plane 5, A (see figure 4): 



\*-^\<Vc*J—^. (34) 
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Consequently, one can always excite both the modes, compensating one detuning by another. For 



example, if A = 45 and v 3> 1, then b 
case equation (|i~2|) yields 



c with |c| <C a 2 , so that \fb ~ |a| — c/(2\a\). In this 



(y + 2J 2 



(35) 



i/ + 25 2 

As far as 5 2 <C f , the maximal value of the increment coefficient Re(A±) is practically the same 
as in the case of strict resonance (|T7|), Energies of both modes have the same order of magnitude 
(as in the strict resonance case), therefore this regime of excitation can be named "symmetrical". 
However, with increase of \S\ the increment of energy growth decreases (by the same law as the 



resonance width r/ c ), and for 5 2 S> v we have Re(A 



±) 



u/(85 2 



1/2 



<C e. 



The second possibility to swing the modes is to choose b < in (14). Then c > a 2 , and 
the eigenvalue A+ is real (whereas A_ is pure imaginary). There are two symmetrical regions of 
swinging in the parameter plane 5, A, which are located between the branches of hyperbolas (see 
figure 4): 



2(1 - 5) 



=- < 7 < 



2(1-5) 



< 7 or 7 < 



2(1 + 5) 
v 

2(1 + 8) ' 



\5\>1, 
\8\<1, 



(36) 
(37) 



where 7 = A — 35. If \5\ < 1, then we may parametrize 7 as 7 = —vt^/4, where |£| > 1. In this 
case, the coefficient c has the following structure: c = 7 4 (1 + 2u/j 2 + C2/7 2 + ...), where C2 ~ 1. 
Using formula (11) and expanding ^[c as \fc = 7 [1 + 1//7 2 + c 2 /(2 7 2 ) - z^ 2 /(2 7 4 ) + ...], we find, 
neglecting corrections of the order of is^ 1 , 



\ 2 + /e 2 = l-5 2 + 4{5£ - l)/e 2 



(38) 



If |^| 3> 1, then we have the same expression for A+ as in the case of a single resonance mode [0,0: 
A + = Ve 2 -^ 2 . This means that the second mode goes out of resonance, if its detuning A satisfies 



the inequality |A| S> v. It is interesting, however, that adjusting two detunings, one can achieve the 
maximal possible value }^ ax ^ ~ e even for 5 ~ e, when the photon generation becomes impossible 
in the single-mode case: this happens for = 2/5, i.e., 7* = —u/(25). 

For \5\ > 1 we write 7 = — \v j{5 + x). Then the resonance exists for — 1 < % < 1. Under the 
condition v 3> 1, the real eigenvalue A+ can be written as (we neglect corrections of the relative 
order of 

A, 



u + 25 2 

Again, one can achieve the maximal possible value X^™ "^ 
maximal eigenvalue decreases as evj(25 2s ). 



(39) 



e, if v > 5 2 , but for v < 5 2 the 
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The resonance width Ta (the distance between two points of intersections between the vertical 
line 5 = const and the hyperbolas limiting the region of resonance in figure 4) equals Ta = vj (5 2 — 1) 
(for \S\ > 1), so it decreases rapidly for \S\ 3> 1. Similar widths Tg (determined in an obvious 
way from the points of intersections between the horizontal line A = const and the hyperbolas) are 

(l r) 

given by the expressions T\ = 1 ± a, where 



a 



(A + 3) 2 + 6z/- V(A-3) 2 + 6 



The sum of these widths does not depend on A, and the smallest of them decreases as Yg Smal1 ^ ps 
3iV(4A 2 ) for |A| > v. 

Simple explicit solutions of the equations of motion can be found in the special case 6 = 1, 
7 = —v/2 (when the generation stops for the single mode, but the increment takes its maximal 
value for coupled two modes). In this case, the eigenvalues of matrix A (|To| ) are as follows, 



A 



+ 



R 



l-- + 0(u 
v 



-2\ 



X_ = iJ, J= - + 1 + 0{u- 2 ). 
With an accuracy up to terms of the order of v~ l , the solutions read as 



x l {t) = xi{0) 

-Pi(0) 
x 3 (0) 



I--] Or (2Rt; t) + - cos 0i 

V ) V 

2\ 2 
- Sr(2RT;t) - - sin 0i 

V V 



+ 



4 M 

p 3 (0) 
12/i 



C{{2Rr;t) -cos0i 
S{(2Rr;t) + sin0i 



(40) 



x 3 (t) = x 3 (0) 



1 - cos 03 + -C7(2Rr;t) 

V V 



+ ^P3(0) 



1 - - ) sin 03 - -S7(2Rr;t) 



+ 



si(o) r 

12/i 

pi(Q) r 

12/i 



C^(2Rr;t) -cos0 3 
S 3 (2Rr;t) +sin0 3 



where the notation is the same as in the preceding section, and 



0fc(r; i) = ku)t — 2 Jr. 



(41) 



The mean energies of each mode depend on time as follows (we neglect corrections of the order 
of is~ 2 , because they are small; moreover, they increase with time slower than the preserved terms 
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of the order of v 1 ): 



4\ 4 
1 - - J cosh(4i?r) + -V'(r) 



+^ [cosh(4i?T) + 1 - 2ip{r)] 



(42) 



4 4 

- + ~Hr) 
v v 



where 



[cosh(4 J Rr) + 1 - 2^(r)] , 



ip(r) = cosh(2i?r) cos(2Jr). 



(43) 



We see that the energy of the first mode is almost the same that it would be in the case of a single 
mode under the condition of the strict resonance, i.e., 9\ cos!i(4t)/2. The energy of the third mode 
is significantly less than the energy of the first mode, if u ^> 1. Therefore this regime of excitation 
can be named "asymmetrical". For r > 1, £-zj£\ ~ 6/za Note, however, that for the cubical cavity 
with v = 50/3, the energy of the third mode is only three times less than that of the first one. It 
is important, nonetheless, that the rates of increase of the energies of each mode are almost twice 
bigger than they were in the case of the strict resonance discussed in the preceding section. 

For the invariant uncertainty products we obtain, with the same accuracy, almost identical 
expressions: 



V 



9 2 
'1,3 



1,3 



s 



1 - - + -V>(r) 



+— [cosh(4 J Rr) + 1 - 2^(r) 



(44) 



We see a drastic difference from the strict resonance case: now the IUP of each mode increases 
(asymptotically) exponentially with time. For r> 1, each mode appears in a highly mixed quantum 
state, with V 1 = V 3 = 6»i6> 3 exp(4i?r)/(2^). 

The photon distributions in the first and the third modes turn out essentially different (and 
different from the distributions in the strict resonance case). We have shown in the preceding 
section, that the form of the PDF (oscillating or not) depends on the ratio EjT>. For the first 
mode, 8\/T>\ ~ u/{26 3 ) for t> 1. Consequently, the asymptotical PDF oscillates for the initial 
vacuum state and for initial thermal states with not very large mean numbers of photons (less 



than, approximately, z^/4). Also, the squeezing coefficient (26) tends to the finite asymptotical 
value sf° = lim^oo (T>i /£\ ) = (typical dependences s\{r) are given in figure 5). This 

resembles (qualitatively) the behavior of the squeezing coefficient in the one-dimensional cavity 
with an infinite number of coupled modes [10]. On the contrary, since T>^/£ 3 — > #3 for r — > 00, the 
final squeezing coefficient of the third mode coincides with the initial value: S3(0) = S3 = 63 > 1, 
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with very small deviations from this value in the interval < r < oo. The photon distribution 
function of the third mode does not oscillate, being close to the "quasigeometric" distribution (j3l|). 

5 Conclusion 

We have studied the problem of photon creation due to the nonstationary (dynamical) Casimir 
effect in a three- or two-dimensional ideal cavity with an oscillating boundary, in the case when the 
spectrum of eigenfrequencies contains two frequencies, whose ratio is close to 3, and the boundary 
oscillates at the frequency close to the double lowest field eigenfrequency. We have calculated the 
mean energy, squeezing coefficient, invariant uncertainty product, and the photon distribution in 
each of two resonance modes, both under the strict and approximate resonance conditions. We 
have shown that, due to a strong coupling between the resonance modes, it is possible to obtain the 
photon generation even for relatively large detunings between the frequencies of the wall and the 
field, compensating one detuning by another. However, the resonance width decreases with increase 
of detuning, as well as the increment of the exponential growth of the energy. We have considered 
both the vacuum and thermal initial states, having demonstrated not only an amplification of 
the number of created photons due to the initial thermal fluctuations, but also such effects as 
the "purity exchange" between the modes and the transformation of the initial smooth photon 
distribution to an oscillating distribution, typical for strongly squeezed states. Our results show 
that choosing a proper set of parameters of realistic three-dimensional cavities, one could facilitate 
observing the nonstationary Casimir effect, which is a challenge for experimentalists. However, a 
correct account of the influence of damping is still an unsolved problem. 
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Figure 1: Mean energies (normalized by their initial values at r = 0) of the first and the third 
modes of the cubic cavity with v = 50/3 in the case of the strict resonance, versus the dimensionless 
"slow time" r (|20|). The order of the curves in the interval of time between 1.4 and 1.6 is as follows 
(from bottom to top): the first mode for $31 = 1/3 (high-temperature initial states); the first mode 
for 9\ = 63 = 1 (vacuum initial state); the third mode for 9\ = #3 = 1; the third mode for #13 = 3. 



Figure 2: The photon distribution function V(n) for the first mode and the initial vacuum state, 
under the condition of the strict resonance, for the "slow time" r = 57r/(2p) and v = 50/3. 



Figure 3: The photon distribution function V(n) for the first mode and the initial thermal state 
with 6\ = 5 (i.e., for (n(0)} = 2 in the first mode), under the condition of the strict resonance, for 
the "slow time" r = 5vr/(2p) and v = 50/3. 



Figure 4: Regions of photon generation in the plane (5, A) of the detuning parameters (tildes 
have been omitted). The region of "symmetrical generation" | A — 4<5| < r\ c is between two dashed 
(almost straight) lines. Two regions of "asymmetrical generation" are confined with hyperbolas 
drawn as continuous curves. 



Figure 5: The squeezing coefficient (|26| ) versus "slow time" r (|2^) for the first mode in the case 
of "asymmetrical generation", for b = 1, 7 = —v/2, and v = 50/3. The lower curve corresponds to 
the case 6\ = 1 (vacuum initial state). The upper curve corresponds to the case 9\ = 5. 
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